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I.  INTRODUCTION  AND  SUMMARY 


As  part  of  the  DNA  stemming  and  containment  program  for  under¬ 
ground  nuclear  testing,  Stanford  Research  Institute  (SRI)  is  conducting 
laboratory  experiments  to  "study  the  residual  stress  fields  around  explo¬ 
sively  formed  cavities.  The  standard  experiment  involves  casting  a 
12-inch  diameter  sphere  of  rock  matching  grout  (2C4)  around  a  lucite- 
encased  sphere  of  high  explosive  (PETN),  placing  this  grout  sphere  in  a 
pressurized  water  tank  to  simulate  overburden  pressure,  and  detonating 
the  PETN.  While  maintaining  overburden  pressure,  the  sphere  is  then 
hydrofractured  from  the  explosively  formed  cavity.  A  detailed 

description  and  discussion  of  the  experimental  results  may  be  found  in 
M 

cizex  ana  riorence'  ’  . 

3 

Systems,  Science  and  Software  (S  )  was  asked  to  numerically 
simulate  these  experiments  in  order  to  increase  understanding  of  the 
laboratory  results,  and  to  validate  our  capability  to  calculate 
containment  related  phenomena.  One  part  of  our  calcul ational  effort, 
the  numerical  simulation  of  the  high  explosive  detonation  and  the 
subsequent  nonlinear  dynamic  processes  which  result  in  the  formation  of 
a  compressive  residual  stress  field  in  the  grout  surrounding  the 
exploded  cavity,  has  been  reported  by  Rimer  and  Lie^.  Since  the 
major  laboratory  data  consist  of  graphs  of  hydrofracture  pressure  vs. 
volume  of  fluid  pumped  into  the  cavity,  the  only  data  available  to  be 
compared  to  the  dynamic  calculations  consisted  of  final  cavity  radius 
and  the  pressure  and  integrated  impulse  at  a  quartz  gauge  emplaced  in 
the  bottom  of  the  water  tank. 

The  reported  calculations  were  in  agreement  with  this  integrated 
impulse  to  within  15  percent  which  is  better  than  the  reproducibility  of 
the  laboratory  record  from  test  to  test.  However,  our  calculated  cavity 
radius  was  approximately  25  percent  higher  than  the  average  measured 
value.  Even  for  this  larger  cavity  radius,  the  calculated  cavity  pres¬ 
sure  before  hydrofracture  was  greater  than  the  average  fracture  initia¬ 
tion  pressure  from  the  laboratory  tests  (a  smaller  calculated  cavity 
radius  would  give  a  still  larger  pressure.)  From  these  results  and  from 
laboratory  experiments  showing  larger  fracture  initiation  pressures  for 
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faster  hydrofractures  we  concluded  that  a  time-related  creep  or  stress 
relaxation  process  must  be  taking  place  in  the  time  between  the  forma¬ 
tion  of  the  residual  stress  field  and  the  hydrofracture  of  the  grout 
spheres.  This  relaxation  process  reduces  the  magnitude  of  these 

stresses  so  that  hydrofracture  can  occur  at  lower  pressures. 

3 

More  recent  unpublished  work  at  S  by  Bill  Proffer  ana  Ed  Peter¬ 
son  has  identified  a  plausible  mechanism  for  this  relaxation  of  the 
residual  stress  field  with  time;  fluid  diffusion  from  saturated  mate¬ 
rials  under  applied  stress  gradients.  A  series  of  calculations  made 
with  the  6ASFL0  code  (a  fluid  diffusion  code)  and  the  SWIS  finite  ele¬ 
ment  code  have  shown  that  the  diffusion  of  water  driven  by  a  residual 
stress  field  in  grout  or  tuff  can  have  a  dramatic  effect  in  reducing  the 
stress  toward  the  ambient.  Within  the  grout  spheres,  the  relaxation 
time  scale  is  measured  in  seconds.  In  the  nuclear  case,  the  calculated 
time  is  of  the  order  of  months.  The  results  suggest  that  the  residual 
stress  fields  being  sampled  in  the  grout  spheres  experiment  are  the 
residuum  left  after  the  above  diffusion  process.  The  nuclear  analog 
would  involve  a  cavity  failure  many  months  after  the  event.  We  strongly 
recommend  that  a  way  be  found  to  speed  up  the  hydrofrac  process  by  an 
order  of  magnitude  or  so  to  test  this  conclusion,  one  possibility  being 
the  use  of  a  gas  pressurization  system  to  hydrofracture  the  exploded 
spheres. 

In  recent  months,  SRI  has  conducted  laboratory  tests  on  configura¬ 
tions  other  than  the  standard  sphere  of  2C4  rock  matching  grout.  Tests 
have  been  made  in  which  the  PETN  sphere  was  placed  inside  a  larger 
cavity  in  the  2C4  grout  sphere  to  simulate  a  decoupled  event.  Tests  are 
also  being  conducted  for  two  different  materials,  a  granite  simulant 
(GS3)  and  a  high  void  content,  low  density  rock  matching  grout  (LD2C4) 
formed  by  adding  a  specified  percentage  of  hollow  silicon  microspheres 
to  the  usual  2C4  mix.  Here,  we  present  the  results  of  our  one¬ 
dimensional,  spherically  symmetric,  dynamic  calculations  of  the  detona¬ 
tion  of  the  PETN  explosive  and  the  formation  of  the  compressive  residual 
stress  fields  around  the  explosive  cavity  for  these  tests. 


Calculations  were  made  for  the  four  basic  test  configurations;  (1) 
the  standard  2C4  sphere,  (2)  the  decoupling  shot,  (3)  the  low  density 
rock  matching  grout  (LD2C4)  and  (4)  the  granite  Simula  it  sphere  (GS3). 
All  of  the  calculations  gave  cavity  radii  significantly  larger  than  the 
measured  values  except  for  the  decoupled  shot  which  gave  reasonable 
agreement.  The  low-density  grout  calculation  showed  the  largest  cavity 
radius,  the  smallest  cavity  pressure,  and  the  smallest  maximum  residual 
hoop  stress.  The  granite  simulant  calculation  gave  the  smallest  cavity, 
but  the  largest  cavity  pressure  and  residual  stress  (cavity  pressure  was 
greater  than  1  Kbar).  In  all  three  coupled  shots,  the  maximum  residual 
hoop  stress  was  approximately  a  factor  of  two  greater  than  the  calcu¬ 
lated  cavity  pressure. 

For  the  decoupled  cavity  shot  the  cavity  pressure  was  about  a 
factor  of  two  greater  than  for  the  coupled  2C4  test.  However,  the  re¬ 
sidual  hoop  stress  for  the  decoupled  cavity  test  was  largest  near  the 
cavity  and  nowhere  exceeded  the  cavity  pressure,  not  an  ideal  situation 
for  containment.  In  spite  of  this,  the  measure d  hydrofracture  pressures 
for  the  uncoupled  tests  were  essentially  the  same  as  for  the  coupled  2C4 
experiments.  This  seems  to  reinforce  the  hypothesis  that  the  residual 
stresses  have  relaxed  with  time.  However,  another  possibility  exists. 
Rimer  and  Friedman^  have  reported  calculations  for  nuclear  events  in 
NTS  Area  12  tuff  in  which  the  results  for  coupled  and  decoupled  events 
are  qualitatively  similar  to  the  results  presented  here.  These  calcula¬ 
tions  also  showed  that  when  the  cavity  pressure  is  allowed  to  decay  to 
zero,  the  peak  residual  hoop  stress  for  the  coupled  and  uncoupled  cases 
are  very  similar  and  only  slightly  smaller  than  for  the  coupled  case 
before  pressure  decay. 

One  disturbing  aspect  of  the  calculations  is  the  larger  cavity 
radii  (27  percent  greater  than  the  measurements  for  the  standard  2C4 
sphere).  A  careful  examination  was  made  of  the  physical  dimensions  of 
the  experiment  including  charge  weight  and  density,  thickness  of  lucite 
shell  around  the  explosive,  location  of  charge  relative  to  the  quartz 
gauge,  and  material  properties  of  the  2C4  grout.  As  a  result,  the 
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uv 6*  A, 


(3) 

calculational  model  described  in  Rimer  ana  Lie'1  '  was  altered;  the 
major  changes  being  a  reduction  in  weight  of  PETN  from  the  nominal  1.375 
gms  to  1.361  gms,  using  slightly  more  lucite,  a  somewhat  different  fail¬ 
ure  envelope  for  low  pressures,  and  a  new  location  for  the  quartz  gauge 
at  a  radius  of  6-5/8  inches  (16.83  cm)  instead  of  at  7  inches  (17.78 
cm).  The  calculated  cavity  radius  decreased  only  about  2  percent  as  a 
result  of  these  changes. 

The  effect  of  moving  the  quartz  gauge  closer  in  was  to  increase 
the  calculated  impulse  by  42  percent,  making  it  approximately  21  percent 
greater  than  the  largest  reported  value  and  31  percent  greater  than  the 
average  of  the  three  reported  values.  Cizek^  estimates  that  the 
gauge  is  located  at  6-5/8  inches  to  within  a  tolerance  of  ^  1/8  to 
1/4-inch.  Based  on  our  calculations  alone,  this  implies  that  the 
possible  error  in  impulse  due  to  locating  the  gauge  alone  may  be  as 
great  as  15  to  31  percent,  irrespective  of  the  reproducibi  1  ity  of  the 
PETN  detonation. 

SRI  has  recently  performed  mechanical  properties  tests  on  the 
three  materials  of  interest,  2C4,  LD2C4,  and  GS3.  Tests  were  conducted 
of  uniaxial  compressive  strength  and  splitting  tensile  strength  under 
both  static  and  dynamic  (strain  rate  of  1.15  or  1.25  sec~^)  loading. 
For  all  three  materials,  the  tests  under  dynamic  loading  conditions 
showed  significantly  higher  compressive  strengths  (31  to  61  percent)  and 
tensile  strengths  (61  to  81  percent). 

5  -1 

Our  calculations  showed  strain  rates  as  large  as  2  x  11  sec 
in  the  sphere  of  2C4  grout,  several  orders  of  magnitude  greater  than 
would  be  expected  in  a  prototype  nuclear  or  HE  event  at  the  same  stress 
levels.  The  laboratory  mechanical  properties  tests  suggest  that  the 
grout  spheres  results  cannot  be  scaled  to  the  large  length  scales  of 
interest  in  the  nuclear  case  because  of  their  large  strain  rates.  Our 
constitutive  strength  models  were  changed  to  take  account  of  the  in¬ 
crease  in  strength  of  the  spheres  under  higher  strain  rate  loading. 
This  required  extrapolating  the  laboratory  strain  rate  data  to  the 
considerably  larger  strain  rates  of  interest.  Calculations  made  using 


this  crude  dynamic  strength  model  gave  cavity  radii  within  5  percent  of 
the  data  for  all  cases  except  the  low-density  grout  (LD2C4)  in  which  the 
calculated  cavity  radius  was  low  by  approximately  10  percent. 

The  larger  strengths  implied  by  the  dynamic  strength  model  re¬ 
sulted  in  a  significantly  narrower  positive  overpressure  pulse  at  the 
quartz  gauge  than  obtained  with  the  original  strength  model  and  thus 
smaller  peak  impulse.  The  peak  pressures  were  basically  the  same. 
Quartz  gauge  data  are  available  only  for  the  standard  2C4  grout  tests 
and  the  low  density  LD2C4  tests.  In  these  cases,  the  measured  durations 
of  positive  overpressure  were  in  far  better  agreement  with  the  calcula¬ 
tions  using  the  dynamic  strength  model  than  the  static  strength  model. 
All  calculated  peak  pressures  were  lower  than  the  data.  The  calcula¬ 
tions  using  the  static  strength  model  gave  impulse  at  the  gauges  larger 
than  the  average  data  by  approximately  10  to  30  percent,  but  for  the 
dynamic  strength  model,  the  impulse  was  smaller  than  the  data  by  30  per¬ 
cent  for  2C4  but  in  good  agreement  for  LD2C4.  Since  the  possible  error 
in  impulse  due  to  locating  the  gauge  only  to  within  _+l/8  to  1 /4-inch  is 
calculated  to  be  15  to  30  percent,  we  feel  that  too  few  measurements 
have  been  made  to  evaluate  our  calculations  in  terms  of  impulse. 

In  summary,  the  calculations  have  demonstrated  that  a  higher 
strength  than  given  by  the  static  laboratory  measurements  may  be  needed 
to  match  the  cavity  radius  data.  Our  crude  model  based  on  the  measured 
dynamic  strengths  gives  a  reasonable  match  to  these  data.  The  higher 
strength  results  in  larger  cavity  pressure  and  higher  peak  y'esidual 
stresses.  These  calculated  values  are  even  less  consistent  with  the 
hydrofracture  pressure  records  even  less  believable  unless  some  time 
dependent  creep  or  stress  relocation  explanation  is  invoked.  This 
provides  further  support  for  our  recommendation  of  a  much  faster 
hydrofracture  procedure. 

The  remainder  of  this  report  consists  of  a  presentation  of  the 
results  of  our  calculations.  In  Section  2  we  discuss  our  constitutive 
models  for  the  three  grout  materials  of  interest  with  particular  empha¬ 
sis  on  the  strength  modeling.  We  present  detailed  results  of  our 
calculations  in  Section  3. 


II.  CONSTITUTIVE  MOOELS  AND  MATERIAL  PROPERTIES  DATA 


In  this  section,  we  discuss  the  constitutive  models  and  material 

properties  used  to  simulate  the  behavior  of  the  different  grouts  used. 

More  complete  discussions  of  the  modeling  for  2C4  grout  and  for  PETN 
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explosive,  lucite  and  water  are  contained  in  Rimer  and  Lie'  .  The 
equation  of  state  for  air  used  for  the  decoupled  calculations  is  from 
Alme^ . 

2.1  Laboratory  Data 


Table  1  gives  the  physical  and  mechanical  properties  data  obtained 
from  SRI ^ for  the  three  materials  of  interest,  the  rock-matching 
grout  (2C4),  the  low-density  rock-matching  grout  (LD2C4)  made  by  combin¬ 
ing  the  basic  2C4  mix  with  hollow  glass  microspheres ,  and  the  granite 
simulant  (GS3).  Additional  data  available  from  Terra  Tek,  Inc.  include 
triaxial  compression  curves^  and  uniaxial  strain  curves^  for 
static  loading  for  2C4  and  GS3  only. 

The  uniaxial  compression  curves  for  2C4  imply  0.7  percent  air 

voids  crushed  up  at  approximately  1.4  Kbars  as  discussed  in  Rimer  and 
(3) 

Lie'  .  Both  uniaxial  and  triaxial  data  indicate  a  maximum  strength 
(stress  difference)  of  0.33  Kb  for  2C4  under  static  loading.  We  will 

discuss  the  dynamic  strain  rate  loading  later  in  this  section.  For  GS3, 

the  uniaxial  curves  show  only  0.5  percent  air  voids  crushed  at  0.9  Kb 
compared  with  3.5  percent  from  Table  1  (calculated  from  water  content 
and  percent  saturated  data).  Triaxial  and  uniaxial  curves  at  0.5  Kb 
confining  pressure  show  stress  differences  of  1.65  Kb  and  1.12  Kb, 

respectively,  but  do  not  indicate  a  maximum  strength  value.  Little  data 
are  available  for  LD2C4. 

These  laboratory  data  were  used  to  construct  equations  of  state 
for  materials  2C4,  LD2C4,  and  GS3.  Reasonable  assumptions  were  made 
where  data  were  lacking. 


2C4  =  Rock -Matching  Grout 

LD2C4  =  Low-Density  Rock-Matching  Grout 

GS3  =  Granite  Simulant 

.  psi  =  6854.76  Pa 

.  bar  =  105  Pa 
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More  recently,  Cizek  and  Florence '  '  have  reported  strength 
measurements  whicu  are  the  average  of  a  set  of  laboratory  measurements 
for  each  of  the  three  materials.  In  some  cases,  these  measurements 
differ  from  the  earlier  numbers  shown  in  Table  1.  They  are  reported  in 
Table  2  for  completeness  but  did  not  influence  the  construction  of  our 
equations  of  state.  Following  is  a  summary  of  the  equations  of  state 
used  for  the  three  materials  of  interest. 

2.2  Rock-Matching  Grout  Model  (2C4) 
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Table  3,  reproduced  from  Rimer  and  Lie'  summarizes  the  mate¬ 
rial  properties  data  used  for  the  early  2C4  grout  calculations  while 
Figure  1  gives  the  load-unload  curves  (porous  crushup)  for  this  mate¬ 
rial.  The  strength  parameters  shown  in  Table  3  are  from  a  parabolic 
variation  of  yield  strength  with  P,  defined  by 


where  P  is  the  mean  stress  and  J'3  the  third  deviatoric  stress  invar¬ 
iant.  These  parameters  give  a  maximum  stress  difference  (Y„  +  V  ) 
of  0.33  Kb  reached  at  F  equal  to  0.4  Kb  and  also  correspond  to  an  uncon¬ 
fined  compressive  strength  of  0.279  Kb  (approximately  4000  psi).  Both 
maximum  strength  (stress  difference)  and  unconfined  strength  are  in 
agreement  with  the  laboratory  data  for  2C4.  However,  the  maximum  pos¬ 
sible  tensile  strength  of  0.263  Kb  (in  a  hydrostatic  state  of  stress)  is 
much  greater  than  given  by  the  laboratory  data  (approximately  0.036  Kb). 
The  model  (no  tensile  failure  without  shear  failure  first)  for  other 
stress  states  does  permit  tensile  failure  in  a  principal  direction  at 
lower  tensile  stress  values  (for  example,  approximately  105  bars  for  a 
state  of  pure  shear). 

A  straight  line  variation  between  yield  strength  and  F  up  to  a 
maximum  strength  of  0.33  Kb  was  used  for  the  present  study.  This  new 
failure  surface  given  by 

Y  -  0.055  +  1.6  P  (kb) 
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TABLE  2.  AVERAGE  STRENGTHS  OF  RMG  2C4 ,  LD2C4 ,  and 

( 8 ) 

GS 3  (Reproduced  from  Cizek  and  Florence  ) 


Material 

Average 

Strain  Rate 
(sec-l) 

Compressive 

Strength 

(psi) 

Tensile 

Strength 

(psi) 

RMG  2C4 

Static 

3990 

580 

RMG  2C4 

0.15 

5330 

900 

LD  2C4 

Static 

3200 

520 

LD  2C4 

0.15 

5000 

780 

GS  3 

Static 

9210 

920 

GS  3 

0.15 

15,900 

1620 
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Table  3.  Summary  of  2C4  Grout 

Material  Proper 

ties 

Data . 

Quantity 

Symbol 

Value 

Units 

Longitudinal  sound 

C 

3.27 

km/sec 

speed 

Shear  wave  speed 

Cs 

1.82 

km/ sec 

Density 

Do 

2.15 

cms/cc 

Zero  pressure  bulk 
modulus 

K 

o 

135.57 

kbars 

Shear  modulus 

G 

71.55 

kbars 

Poisson's  ratio 

a 

0.275 

(  A 

150 

kbars 

3 

335 

kbars 

i 

Coefficients  for 

5 

5 

Tiilotscr.  Equation  < 

'  a2 

0.5 

of  State 

)  ’D 

0.633 

6.0x10^ 

3 . 5x10  *■!* 
l.SxlO1- 

1 

I  eo 
e<? 
esl 
v  es 

ergs/gm 

ergs/gm 

ergs/gm 

| 

f 

,  Y 

0.21 

kbars 

Strength  Parameters  \ 

^  j 

1  vc 
(  ?m 

0.12 

0.40 

kbars 

kbars 

1 

em 

'  m 

2x10-0 

ergs/gn 

Air-filled  Porosity 

on 

0  .007 

Elastic  Pressure 

0.1 

kbars 

Crush  Pressure 

*pe 

c 

1.4 

kbars 

Figure  1.  Load-unload  curves  for  2C4  grout 
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is  consistent  with  both  the  laboratory  values  of  unconfined  compressive 
strength  and  tensile  strength  for  2C4  grout.  Figure  2  shows  a  compar¬ 
ison  between  the  old  and  new  failure  surfaces. 

2.3  Granite  Simulant  Model  (GS3) 

Whenever  possible,  laboratory  data  were  used  in  developing  the 
constitutive  model  for  the  granite  simulant.  Thus,  the  density  and 
elastic  parameters  for  GS3  come  directly  from  Table  1  as  does  the  calcu¬ 
lated  air  voids  of  3.5  percent.  The  uniaxial  strain  load-unload 

curve^  showed  only  0.5  percent  air  voids  crushed  up  at  a  mean  stress 
of  0.9  Kb.  Assuming  that  the  calculated  air  voids  indicated  a  consider¬ 
ably  higher  crush  pressure,  we  used  the  0.5  percent  air  void  at  0.9  Kb 
as  one  point  on  the  loading  curve  shown  in  Figure  3.  The  loading  curve 
above  the  crush  pressure  was  a  straight  line  with  bulk  modulus  of  335  Kb. 

In  a  recent  conversation,  0.  C.  Cizek  of  SRI  has  brought  into 

question  the  saturation  data  for  GS3  given  in  Table  1.  The  experimental 
procedure  is  to  maintain  the  GS3  sphere  in  a  water  bath  at  all  times. 
Therefore,  Cizek  believes  that  the  granite  simulant  is  fully  saturated, 
and  that  the  air-filled  voids  are  negligible.  This  informaton  came  too 
late  to  be  included  in  our  constitutive  model  and  calculations.  However, 

we  believe  that  since  GS3  was  not  mixed  in  a  vacuum,  air  bubbles  were 

entrained  which  could  result  in  a  few  percent  air-voids.  The  water  bath 
could  not  remove  such  bubbles. 

If  the  simulant  were  truly  100  percent  saturated,  we  would  antici¬ 
pate  a  slightly  larger  cavity  radius  than  calculated  with  our  model. 
However,  the  measured  cavity  radius  was  actually  smaller  than  calcu¬ 
lated.  The  absence  of  air-filled  voids  would  increase  the  pressure  and 
impulse  at  the  quartz  gauge  greatly.  No  data  are  available  for  compar¬ 
ison. 

The  measured  unconfined  compressive  strength,  tensile  strength  and 
stress  difference  at  0.5  Kb  confining  pressure  under  triaxial  loading 
were  used  to  determine  the  coefficients  of  a  parabolic  failure  surface 
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Figure  2.  Old  and  new  failure  surfaces  for  2C4  grout. 


Figure  3.  Specific  volume  load-unload  curves 
for  granite  simulant  (GS3) . 
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Figure  4.  Failure  surface  for  granite  simulant  GS3 


(vs  P)  for  static  loading  conditions.  Since  the  uniaxial  loading  data 
did  not  show  a  clear  maximum  stress  difference,  and  since  no  other  data 
were  available,  we  extrapolated  the  parabolic  form  of  the  failure  sur¬ 
face  to  obtain  a  maximum  stress  difference  of  1.795  Kb,  only  9  percent 
greater  than  the  given  triaxial  measurement.  Figure  4  shows  our  failure 
surface  for  GS3.  Note  that  the  uniaxial  data  (at  0.5  Kb  confining  pres¬ 
sure)  is  below  the  failure  surface. 

2.4  Low  Density  Rock  Matching  Grout  Model  (LD2C4) 

The  low  density  rock  matching  grout  (LD2C4)  was  made  by  blending 
the  standard  2C4  mix  with  a  specified  amount  of  Q-CEL  microspheres  which 
encapsulate  air  in  hollow  spherical  shells  derived  from  inorganic  sili¬ 
cate.  Our  first  model  utilized  the  sparse  preliminary  data  available  at 
that  time  (13.4  percent  air  voids  and  a  density  of  1.92  gns/cm  )  and 
filled  in  the  gaps  by  using  the  material  properties  developed  for  2C4. 
Thus,  the  crush  pressure  P  ,  the  elastic  pressure  P  and  the  general 
shape  of  the  load-unload  curves  shown  in  Figure  5  are  similar  to  Figure 
1.  The  straight  line  failure  surface  of  Figure  2  was  used  for  the  first 
L02C4  grout  spheres  calculation. 

Our  present  model  incorporates  the  later  data  shown  in  Table  1. 
The  major  changes  from  the  preliminary  model  are  in  the  failure  surface, 
given  by 

Y  =  0.0525  +  1.52  P-  (Kb) 

with  a  specified  maximum  strength  of  0.26  Kb,  which  is  approximately  20 
percent  smaller  than  before,  and  in  the  load-unload  curve.  We  felt  that 
the  crush  pressure  should  be  considerably  larger  because  of  the  large 
amount  of  air-filled  voids.  Retaining  the  same  general  shape  for  the 
load  curve  gave  us  the  present  curve  shown  in  Figure  5. 
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2.5  Dynamic  Strength  Models 

Table  1  gives  the  values  of  unconfined  compressive  strength  and 
tensile  strength  for  the  three  materials  of  interest  for  both  static  and 

dynamic  loading.  (The  average  strength  values  shown  in  Table  2  came  too 

late  to  be  included  in  our  models.)  For  our  purposes,  we  assumed  that 
the  "static"  loading  was  represented  by  an  average  strain  rate  of  0.001 
sec-''’.  Thus,  the  "dynamic"  loadings  correspond  to  approximately  two 
orders  of  magnitude  greater  strain  rates.  At  these  higher  strain  rates, 
the  unconfined  compressive  strength  was  increased  by  30  to  60  percent 
(depending  on  the  particular  material)  while  the  tensile  strength  was 
increased  by  60  to  80  percent.  Since  no  laboratory  data  were  available 

on  increased  shear  strength  at  confining  pressure,  we  decided  to  assume 

that  maximum  strength  (stress  difference)  was  increased  by  the  same  per¬ 
centage  as  the  unconfined  strength. 

The  classical  theory  for  metals  assumes  that  strength  increases 

linearly  with  the  logarithm  of  strain  rate.  Our  grout  spheres  calcula- 

5  -1 

tions  have  shown  strain  rates  of  the  order  of  10  sec  ,  eight 
orders  of  magnitude  greater  than  the  assumed  static  loading  data  and 
approximately  6  orders  of  magnitude  greater  than  the  dynamic  loading 
data  of  Table  1.  Thus  the  strain  rate  model  must  extrapolate  the  exist¬ 
ing  data  over  many  orders  of  magnitude  and,  therefore,  could  be  quite 
inaccurate.  The  difference  between  using  the  data  of  Table  1  or  Table  2 
when  extrapolated  to  higher  strain  rates,  could  be  significant.  How¬ 
ever,  we  hoped  that  our  model  would  show  us  whether  or  not  strain  rate 
effects  could  bring  the  calculated  cavity  radii  into  better  agreement 
with  the  measurements . 

Our  model  assumes  that  at  all  strain  rates  the  functional  form  of 
the  failure  surface  (straight  line  for  2C4  and  LD2C4  and  parabolic  for 
GS3)  does  not  vary.  For  a  given  time  and  Lagrangian  position,  the  fail¬ 
ure  surface  (stress  difference  as  a  function  of  P)  is  defined  by  the 
radial  strain  rate  in  the  computational  cell.  Deviatoric  stresses  at 
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this  time  and  position  are  limited  to  the  present  failure  surface  as 
described  in  Rimer  and  Lie  .  At  the  next  instant  of  time,  the  failure 
surface  for  the  cell  will  change  if  the  strain  rate  varies. 

The  constants  defining  the  functional  form  are  themselves  defined 
as  functions  of  strain  rate.  For  example,  we  have  already  defined  a 
simple  failure  surface  for  2C4  grout  under  static  loading  (0.001 
sec-*)  using  three  constants,  a  stress  difference  at  zero  P,  a  slope 

of  the  straight  line  segment,  and  a  maximum  stress  difference.  We  use 
the  data  from  Table  1  for  a  strain  rate  of  0.15  sec-*  to  define  a 
similar  failure  surface  for  that  strain  rate  and  relate  the  three  con¬ 
stants  (functions  of  strain  rate)  by  assuming  linear  variations  of  these 

constants  with  logarithm  of  the  strain  rate.  For  the  granite  simulant 
6S3,  the  procedure  is  similar  except  that  a  parabolic  form  is  assumed. 
Table  4  defines  the  strain  rate  dependence  for  the  three  materials  of 
interest. 

III.  CALCULATIONAL  RESULTS 
3.1  Description  of  Calculations 

Table  5  describes  and  summarizes  the  results  of  the  calculations 
made  using  the  spherically  symmetric  SKIPPER  code.  In  this  table, 
static  failure  refers  to  the  standard  failure  model  based  on  data  from 
laboratory  tests  under  static  loading  conditions  while  dynamic  failure 
refers  to  the  model  presented  in  section  2.5  which  used  the  laboratory 

data  under  dynamic  loading  conditions  presented  in  Table  1.  Calcula¬ 
tions  1  through  4  were  made  using  the  model  for  the  explosive  lucite 

3 

sphere  described  in  Rimer  and  Lie  ,  i.e.,  exactly  0.375  gm  of  PETN, 
outer  radius  of  0.447  cm,  surrounded  by  lucite  out  to  a  radius  of  0.476 
cm.  All  calculations  were  for  a  grout  sphere  15.24  cm  in  radius,  sur¬ 
rounded  by  water  at  a  pressure  of  1000  psi  (approximately  69  bars).  For 
the  first  four  calculations,  the  quartz  gauge  was  modeled  as  being  7 
inches  from  the  center  of  the  sphere.  For  all  other  calculations,  the 
gauge  was  at  6-5/8  inches. 
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TABLE  4.  DEPENDENCE  OF  FAILURE  SURFACES  ON  STRAIN  RATE 


2C4 

Y  =  0.055  f1  +  1.60  f2  F  (Kb) 
Ymax  -  0.33  f3  (kb) 

f,  =  1  +  0.32  7 

f2  -  1  -  0.0287  7 
f3  =  1  +  0.155  7 

where 

7  =  log-jg  (strain  fate/0.001) 


LD2C4 

V  =  0.0525  f4  +  1.52  P  (Kb) 

Ymax  *  °*26  f5  <Kb> 
f4  =  1 .  +  0.245  c 
f5  -  1.  ♦  0.237  7 


GS3 


Y  a 

Y 

max 


0.115  fg 
=*  0.115 


1  _  r  PI  — 

+  T38f?  P  [2  “  1 .88f“7]  (Kb^  P  < 
f6  +  1.68  fg  (Kb)  7  >  1.88  (Kb) 


1.88  (Kb) 


fg  -  1  +  0.317  e 

f7  «  1  +  0.483  7 
fg  -  1  +  0.439  7 
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As  described  in  Section  2.2.,  the  material  properties  used  in 
Calculation  1  differed  from  those  of  Reference  3  only  in  the  shape  of 
the  failure  surface  for  2C4  grout.  The  results  of  this  calculation  were 
essentially  the  same  as  reported  earlier  (peak  reflected  impulse  in¬ 
creased  from  735  to  766  bars-usec).  Calculation  2  used  the  same  mate¬ 
rial  properties  but  modeled  air  between  the  lucite  shell  and  the  2C4 
grout  at  a  radius  of  0.95  cm.  For  this  decoupled  test,  we  initialized 
the  stress  distribution  in  the  grout  using  the  lame  solution  for  a 
hollow  elastic  sphere  pressurized  from  the  outside.  Calculation  3  used 
the  material  properties  for  the  granite  simulant  GS3  discussed  in 
Section  2.3  while  Calculation  4  incorporated  the  preliminary  model  for 
L02C4  grout  from  Section  2.4. 

For  Calculation  5,  the  mass  of  PETN  explosive  was  reduced  slightly 
to  0.360  gm  and  the  lucite  sphere  was  given  an  outer  radius  of  0.482  cm 
reflecting  more  current  information  about  the  SRI  tests.  These  changes 
only  insignificantly  altered  the  calculational  results  in  the  grout  when 
compared  with  Calculation  1.  However,  another  change,  moving  the  simu¬ 
lated  location  of  the  quartz  gauge  in  to  6-5/8  inches  (16.8  cm),  gave  a 
large  increase  (42  percent)  in  the  calculated  reflected  impulse.  (This 
result  will  be  discussed  in  some  detail  later  in  this  report.)  All  sub¬ 
sequent  calculations  used  this  gauge  location  and  the  charge  configura¬ 
tion  of  Calculation  5.  Calculation  5  used  more  recent  laboratory  data 
to  better  define  the  constitutive  model  for  LD2C4  grout.  These  resulted 
in  a  lower  shear  strength  and  less  easily  crushable  air-filled  voids. 

Calculations  2,  3,  5  and  6  include  our  best  estimates  at  present 
of  the  constitutive  models  and  material  properties  for  the  tests  of 
interest  based  on  the  standard  laboratory  materials  tests.  Because  of 
the  lack  of  agreement  with  some  of  the  SRI  test  results  the  constitutive 
model  described  in  Section  2.5  was  developed  to  simulate  the  effect  of 
dynamic  (i.e.,  strain  rate  dependent),  loading  on  the  shear  failure  of 
the  grouts  tested.  Calculations  7  through  10  of  Table  5  were  made  using 
this  "dynamic  failure"  model. 
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3.2  Comparisons  Between  Calculations  and  Cavity  Radius  and  Quartz 

Gauge  Data  " 

Unfortunately  very  little  experimental  data  exist  which  is  di¬ 
rectly  comparable  to  the  calculations.  We  can  compare  calculated  and 
measured  cavity  radii  for  all  four  of  the  experimental  configurations  of 
interest,  cOupled  and  decoupled  2C4  grout,  granite  simulant  GS3  and  low 
density  grout,  LD2C4.  However,  quartz  gauge  data  are  only  available  at 
the  present  time  for  the  2C4  and  L02C4  grout  configurations. 

Table  5  summarizes  many  of  the  calculational  results  and  compares 
them  with  the  measurements.  In  general,  for  the  static  failure  model, 
(Calculations  1  through  6)  calculated  cavity  radii  are  considerably 
greater  than  the  measured  values,  the  only  exception  being  the  decoupled 
2C4  cavity  where  the  calculated  radius  is  within  the  range  of  the  data. 
Using  our  strain  rate  dependent  failure  surface  to  increase  the  shear 
strength  (Calculations  7  through  10)  in  all  cases  resulted  in  a  lower 
calculated  cavity  radii  than  for  the  failure  surfaces  based  on  static 
laboratory  tests.  The  three  coupled  calculational  configurations  gave 
decreases  of  approximately  20  percent  in  cavity  radius  while  the  de¬ 
coupled  2C4  cavity  (where  shear  failure  is  relatively  insignificant) 
decreased  in  radius  less  than  8  percent.  For  both  the  coupled  2C4  and 
the  GS3  configurations,  the  resulting  cavity  radii  were  within  5  percent 
of  the  experimental  measurements.  However,  the  LD2C4  grout  calculations 
which  showed  the  largest  decrease  in  cavity  radii  due  to  dynamic  effcts, 
gave  significantly  smaller  radius  than  the  data.  This  may  be  due  to  the 
extrapolation  of  strength  data  based  on  Table  1  rather  than  the  later 
data  of  Table  2.  Taken  as  a  whole,  the  results  using  our  crude  dynamic 
failure  model,  which  extrapolates  the  laboratory  data  many  orders  of 
magnitude  in  strain  rate,  indicate  the  importance  of  strain  rate  depend¬ 
ent  effects  in  these  small  scale  laboratory  experiments.  Our  guess  is 
that  these  effects  would  be  insignificant  for  the  prototype  nuclear 
tests  in  which  strain  rates  may  be  three  orders  of  magnitude  lower. 

Figure  6  shows  the  pressure-time  record  at  the  quartz  gauge  in  the 
floor  of  the  water  tank  for  2C4  grout  Test  157  together  with  the  simu¬ 
lated  records  for  Calculations  1,  5,  and  7.  Test  157,  the  most  recent 
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of  the  three  available  quartz  gauge  records,  gave  the  smallest  inte¬ 
grated  impulse  of  751  bars-usec  (the  maximum  was  896  bars-vsec).  The 
gauge  is  nominally  located  6-5/8  inches  from  the  center  of  the  explosive 
but  only  to  within  +  1/8  to  1/4  inch.  Calculation  I  placed  this  gauge 
at  7  inches  while  Calculations  5  and  7  located  it  at  6-5/8  inches.  All 
three  calculations  show  an  earlier  signal  arrival  due  to  numerical  dif¬ 
fusion.  A  more  valid  comparison  is  between  the  times  of  arrivals  of  the 
peak  (reflected)  pressures.  Calculation  1  has  the  peak  arriving  too 
late;  Calculation  5  has  the  same  arrival  as  the  data;  and  Calculation  7 
(the  dynamic  failure  model)  has  a  slightly  early  arrival.  All  calcula¬ 
tions  for  2C4  grout  show  a  lower  peak  pressure  than  the  gauge  record. 
However,  the  shapes  of  the  calculated  pulses  and  the  impulse  derived 
from  them  vary  greatly.  Calculation  1  has  the  same  shape  as  the  charac¬ 
teristic  gauge  record,  i.e.,  a  rise  to  the  peak  followed  by  a  fall  to  a 
plateau  and  a  further  fall-off.  Calculation  5,  for  the  closer-in  gauge 
location  and  the  slightly  different  charge  configuration,  shows  a  rise 
to  a  second  peak  following  the  fall-off  which  appears  to  be  due  to  a 
small  second  pulse  (present  in  Calculation  1  as  well)  which  is  in  phase 
with  the  reflected  main  pulse  at  the  gauge  for  this  particular  gauge 
location.  This  small  pulse  appears  to  be  a  reflection  from  the  lucite 
PETN  interface  which  has  rebounded  at  the  origin  and  is  propagating 
radially  outward.  Neither  a  plateau  nor  second  peak  is  visible  in  the 
pulse  from  Calculation  7  for  the  dynamic  failure  model.  The  pressure 
record  from  Calculation  7  would  give  an  impulse  in  excellent  agreement 
with  the  data  if  a  plateau  such  as  in  Calculation  1  was  present.  Table 
5  shows  comparisons  between  the  simulated  and  experimentally  derived 
impulses.  Calculation  1  gives  excellent  agreement  (it  has  a  wider  pulse 
but  lower  peak)  with  the  data,  while  Calculation  5  gives  significantly 
larger  impulse  due  to  the  second  peak. 

Figure  7  compares  the  quartz  gauge  records  from  tests  201  and  202 
with  three  calculations  for  L02C4  grout,  the  preliminary  model  (Calcula¬ 
tion  4)  with  gauge  location  at  7  inches,  the  later  model  with  static 
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Figure  7.  Comparison  between  quartz  gauge  records  (Tests  201  and  202) 
and  calculations  for  LD2C4  group. 


failure  data  (Calculation  6)  and  with  dynamic  failure  data  (Calculation 
10)  both  with  the  gauge  location  at  6-5/8  inches.  As  was  true  for  2C4 
grout  comparisons,  arrivals  are  earlier  than  the  data.  Both  Calcula¬ 
tions  4  and  6  (the  static  failure  model)  show  much  larger  pulse  widths 
than  the  gauge  records.  However,  the  dynamic  failure  model  gave  a  pulse 
width  in  excellent  agreement  with  the  data.  Again,  the  dynamic  model 
gave  a  pulse  shaped  differently  than  the  data.  Tests  201  and  202  (the 

early  records  sent  to  us)  used  spheres  of  LD2C4  grout  that  had  higher 

(9) 

than  normal  densities.  Cizek  describes  the  standard  LD2C4  mix 
(density  1.9)  as  giving  considerably  lower  peak  pressures  (17  bars)  and 
maximum  impulses  (200  bars-usec)  than  these  tests  and  having  a  shape 
similar  to  the  2C4  records.  The  pulse  of  Calculation  10  (the  strain- 
rate  dependent  model)  is  in  excellent  quantitative  agreement  with  the 
results  from  the  tests  using  the  standard  LD2C4  mix. 

Figure  8  shows  the  calculated  pressure  vs  time  for  the  granite 
simulant  6S3  for  the  static  model  (gauge  at  7  inches)  and  the  dynamic 
failure  model  (gauge  at  6-5/8  inches).  Figure  9  shows  similar  curves 
for  the  decoupled  2C4  calculations.  In  both  figures,  the  only  notice¬ 
able  differences  between  static  and  dynamic  failure  models  are  the 
earlier  arrivals  due  to  change  in  gauge  location  and  the  narrower  pulses 
for  the  dynamic  models.  No  data  are  available  for  comparison. 

3.3 _ Discussion  of  Results 


The  major  purpose  of  the  grout  spheres  tests  was  to  evaluate  the 
importance  of  the  explosively  produced  compressive  residual  stress 
fields  on  containment  of  the  cavity  gases.  We  believe  that  due  to  a 
time  dependent  stress  relaxation  mechanism  the  residual  stresses  of 
interest  are  no  longer  present  in  the  laboratory  experiments  at  the  time 
of  hydrofracture.  However,  we  hope  that  our  calculations,  together  with 
the  existing  laboratory  data,  may  prove  useful  in  comparing  containment 
prospects  for  the  four  laboratory  configurations,  in  order  to  extrap¬ 
olate  the  results  to  nuclear  test  media. 
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Pressure  (bats) 


Figure  8 .  Simulated  quartz  gauge  records  for 
granite  simulant  GS 3  tests. 


I*  res 


Figure  10  gives  the  calculated  peak  stresses  (relative  to  the  69 
bar  overburden)  vs  range  from  all  of  the  calculations  for  the  static 
model  (Calculations  1  and  5  overlay)  and  from  Calculation  10  for  L02C4. 
The  plots  for  static  and  dynamic  loading  for  all  configurations  were 
very  similar  down  to  0.5  Kb  or  less  when  strength  effects  become  signif¬ 
icant.  Since  the  magnitude  of  residual  stress  fields  has  been  shown 
(see  Rimer^^)  to  be  very  dependent  on  yield  strength,  peak  stress 
alone  is  not  a  good  indication  of  containment.  The  situation  is  further 
complicated  by  the  effects  of  other  material  properties  on  cavity  size, 
cavity  pressure,  and  stress  levels. 

Let  us  compare  the  results  of  Table  5  and  Figure  10  first  for 
LD2C4  grout.  Calculation  4  (the  preliminary  model)  and  Calculation  6 
(the  later  model)  show  very  similar  peak  stresses  down  to  approximately 
10  Kbar  (the  higher  stress  levels  are  dominated  by  air-filled  voids 
which  are  the  same  for  both).  However,  in  the  range  between  1  and  10 
Kb,  the  results  differ  because  of  the  different  crush  pressures  used. 
Calculation  6  shows  a  noticeable  increase  then  decrease  in  slope  in  the 
peak  stress  curve  at  a  stress  level  of  5  Kb  (the  crush  pressure  is  4  Kb) 
while  Calculation  4,  for  a  crush  pressure  of  1.4  Kb  (also  a  20  percent 
lower  strength),  shows  this  phenomena  at  a  lower  stress  level.  The 
result  is  considerab ly  more  air-voids  crushed  up  for  Calculation  4  and 
lower  peak  stress,  impulse,  etc.  Even  though  more  voids  are  crushed, 
the  cavity  is  slightly  smaller  (and  cavity  pressure  is  higher)  due  to 
the  higher  strength  used  in  Calculation  4  so  that  the  ratio  of  peak 
residual  hoop  stress  to  cavity  pressure,  assumed  here  to  be  a  rough 
containment  guide,  is  lower.  Calculation  10,  identical  to  Calculation  6 
except  for  the  stronger  failure  surface,  shows  a  factor  of  2  higher 
residual  stresses  but  also  a  greater  than  factor  of  2  increase  in  cavity 
pressure.  In  fact,  the  trend  for  all  tamped  configurations  is  that  the 
dynamic  model  results  in  lower  ratios  of  peak  residual  hoop  stress  to 
cavity  pressures . 

Primarily  due  to  the  lover  air-void  content,  the  2C4  calculations 
show  considerably  larger  peak  stresses  at  all  ranges  and  smaller  cavi¬ 
ties  than  for  L02C4.  Approximately  the  same  ratios  of  peak  residual 
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hoop  stress  to  cavity  pressure  are  seen.  The  granite  simulant  G S3 

sphere  has  a  volume  of  air-filled  voids  intermediate  between  2C4  and 
L02C4  so  that  the  peak  stress  vs  range  curve  lies  in  between  the  other 
two  conf i gurations  over  the  stress  regime  where  the  effects  of  voids 
dominate.  Once  again,  changes  in  slope  are  seen  at  peak  stresses  around 
the  crush  pressure  of  3  Kb,  but  at  stress  levels  below  0.3  Kb  the  curve 
is  almost  identical  to  the  2C4  curve.  The  high  strength  of  GS3  results 
in  very  small  cavities,  very  large  cavity  pressures,  and  very  high 

residual  stresses  when  compared  with  the  other  configurations.  However, 
the  ratio  of  peak  residual  hoop  stress  to  cavity  pressure  is  smaller, 
which  may  indicate  worse  containment  prospects. 

Calculations  2  and  8,  for  the  decoupled  configuration,  give  some¬ 
what  dicferent  results.  !n  Figure  10,  peak  stresses  are  much  higher  in 
close  due  to  the  impact  of  the  shocked  air  on  the  walls  of  the  decoupled 
cavity  and  attenuation  is  much  steeper.  The  cavity  expands  only  slight¬ 
ly  from  its  initial  radius  of  0.95  cm.  As  a  result,  considerably  less 
plastic  yielding  occurs  around  it. 

Figures  11  and  12  show  the  residual  stress  fields  calculated  for 
the  decoupled  2C4  spheres  using  the  static  and  dynamic  failure  models, 
respectively.  Figures  13  and  14  show  similar  plots  calculated  for  the 

coupled  2C4  spheres  while  Figures  15,  16,  17  and  18  show  plots  for  the 

coupled  L02C4  and  GS3  spheres.  The  times  at  which  the  stress  fields  are 
monitored  are  included  since  elastic  rebounds  due  to  the  simulated  tank 
wall  are  present  which  modify  these  fields  out  to  much  later  times. 
(These  rebounds  may  change  the  peak  residual  stresses  from  the  values 
shown  by  _+10  percent.)  In  all  cases  studied,  the  peak  residual  stresses 
are  both  greater  for  the  dynamic  failure  model  than  for  the  static  mrdel 
and  are  located  at  smaller  radial  distances.  The  coupled  grout  spheres 
show  peak  residual  hoop  stress  significantly  greater  than  the  peak 
radial  stress.  However,  the  decoupled  calculations,  due  to  the  smaller 
cavity  expansion  and  less  plastic  yielding,  give  residual  stresses  that 
appear  to  be  more  like  the  elastic  solution,  i.e.,  both  maximum  hoop 
stress  and  radial  stress  located  at  or  near  cavity  boundary  and  radial 
stress  greater  than  hoop  stress  over  the  ranges  of 
interest. 
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Radial  Distance  icm) 


Figure  12.  Residual  stress  fields  at  100  ysec  for 
decoupled  2C4  calculation  8  (dynamic 
failure  model) . 
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Figure  13.  Residual  stress  fields  at  84  usee  for  2C4 
calculation  1  (static  failure  model) . 
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Figure  14.  Residual  stress  fields  at  79  usee  for  2C4 
calculation  7  (dynamic  failure  model) . 
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Figure  15.  Residual  stress  fields  at  124  usee  for 

LD2C4  calculation  6  (static  failure  model) 
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Residual  stress  fields  at  95  usee  for  GS3 
calculation  3  (static  failure  model) . 


(hens) 


4000 
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Figure  18.  Residual  stress  fields  at  74  usee  for  GS3 
calculation  9  (dynamic  failure  model) . 
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